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ABSTRACT 

We study the relation between the 300-700 Hz upper kHz quasi-periodic oscillation (QPO) and the 
401 Hz coherent pulsations across all outbursts of the accreting millisecond X-ray pulsar SAX J1808.4- 
3658 observed with the Rossi X-ray Timing Explorer. We find that the pulse amplitude systematically 
changes by a factor of ~ 2 when the upper kHz QPO frequency passes through 401 Hz: it halves when 
the QPO moves to above the spin frequency and doubles again on the way back. This establishes for 
the first time the existence of a direct effect of kHz QPOs on the millisecond pulsations and provides 
a new clue to the origin of the upper kHz QPO. We discuss several scenarios and conclude that while 
more complex explanations can not formally be excluded, our result strongly suggests that the QPO 
is produced by azimuthal motion at the inner edge of the accretion disk, most likely orbital motion. 
Depending on whether this azimuthal motion is faster or slower than the spin, the plasma then 
interacts differently with the neutron-star magnetic field. The most straightforward interpretation 
involves magnetospheric centrifugal inhibition of the accretion flow that sets in when the upper kHz 
QPO becomes slower than the spin. 

Subject headings: pulsars: general - stars: neutron - X-rays: binaries - individual (SAX JI808.4-3658) 


1. INTRODUCTION 

The low-mass X-ray binary SAX JI808.4-3658 
(SAX J1808) is a 401 Hz accreting millisecond X-ray 
pulsar (AMXP; iWi.ina nds fc van der Kli^|1998|) with 
a 2 hr orbital period (Chakrabarty & Morgan’1998). 


The system is an X-ray transient and has shown hve 


X-ray outbursts between 1 998 and 2011 (Hartman et al. 


20081 iPatruno et al. [20121) that were observed with the 

k 


Rossi X-ray Timing Explorer (RATE). The pulsations of 
SAX J1808 are near-s i nusoidal with amplitudes o f 2-7% 
(Hartman et al.||2008 2009| Patruno et al.||2012|. They 
are produced by thermal hotspots, emission regions near 
the stellar magnetic poles, under going periodic aspect 
varia t ions due to the stellar spin (Davidson & Ostriker 
1973 Ghosh & Lamb 1978). The h otspots, which may 


cover a large traction of the surface (Lamb et al. 
are heated by the impact of plasma channeled 


2009), 
ly the 


magnetic held towards the poles. The pulsations there¬ 
fore offer a direct probe of the inner accretion how and 
specihcally of the magnetic threading of the accretion 
disk. 


SAX J1808 was th e hrst AMXP discovered (Wijnands 
& van der Klis 1998), and also the hrst pulsar t o show tEe 
twin kHz quasi-p'eriodic oscillations (QPOs; Wijnands 


et al. 2003) now known to be ubiq uitous among accre ting 
low-inagnetic held neutron stars (van der Klis 2006). In 
SAX J1808 usually only the higher-trequencyof the twin 


et al. 2005|), at frequencies up to 700 Hz 


peaks is observed (jWiinands et ahl 2003 

70 ^ 


van Straaten 
and it is this 


upper kHz QPO we report on here. The lower kHz Q PO 
is seen only on rare occasions ( [Wijnands et al.|[2003 ). 
Some kHz QPO models explain the upper kHz QP 


QPO models explain the upper kHz QPO 
with orbital motion of short-lived inhom ogeneities at 
some preferred radius in the accretion how (|van der Klis 
et al. 1996[ [Strohmayer et al. 1996 iVliller et al. 199^ 


Stella & Vietri 19991 |Kluzniak et al. 20041 Alpar 

& Psaltis 

20081), others invoke alternative mechanisms 1 

Lai 1998 

Kato 

2004 

Zhang 

2004 

Bachetti et al. 

2010 

), but tor 


remained elusive. 

Although both the plasma channeling responsible for 
the X-ray pulsations and the mechanism producing the 
kHz QPOs are believed to probe the inner accretion how 
of AMXPs, no direct relation between pulsations and 
kHz QPOs has been observed so far. In this Letter we 
show that such a direct relation does in fact exists in 
SAX J1808. This provides a new clue to the origin of 
the upper kHz QPO and the nature of the inner accretion 
how towards this AMXP. 

2. DATA REDUCTION 

We consider all pointed observations of SAX J1 808 
with RATE using the Proportional Counter Array (Ja- 


p.m.bult@uva.nl 


hoda et al. 2006[). Each observation consists of one or 
multiple continuous 3 ks exposures as set by the data 
gaps associated with the 95 minute satellite orbit. We 
use the 16 s time-resolution Standard-2 data to create 
2-16 ke V light curves, which we n ormalize to the Crab 
(see, e.g, van Straaten et al. 2003 for details). 

Eor the stochastic timing analysis we use all 
GoodXenon and Event data in the 2-2 0 keV energy 
range . We compute Leahy normalized (Leahy et al. 
1983|) power spectra, using 256 s data segments binned 
to a 1/8192 s (~ 122 ps) time resolution, giving a fre¬ 
quency resolution of 1/256 Hz and a Nyquist frequency 
of 4096 Hz. No background subtraction or dead-time 
correction was done before a power spectrum was calcu¬ 
lated. The resulting power spectra were averaged over 
intervals of 1-15 ks (depending on the number of active 
PCUs, which was higher for observations early in the 
RATE mission) as required to signihcantly detect the 
upper kHz QPO. 
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We subtracted a modeled Poisson noise power spec¬ 


trum (Zhang et al. ] j)95) from the averaged p ower spectra 


using the method ot Klein-Wolt et al. (20041. Accounting 


for the background emission, we then renormalized the 


(van der Klis 

1995 

). 

We htted t 

re power spectra with a set ot Lorentzian 

profiles (|Belloni et al. 

2002), with each profile 


Ti), a tunction ot Fourier frequency v, de¬ 
scribed by three parameters: the centroid frequency z/o,i, 
the full-width-at-half-maximum Wi, and the fractional 
rms 


Ti = 


Li(y)dv. 


( 1 ) 


The subscript i refers to the specific power spectral com¬ 
ponent; in particular, represents the upper kHz QPO. 

When its centroid frequency is high, Lu has good co¬ 
herence and appears as a narrow peak in the power spec¬ 
trum; as the frequency decreases, so does the coherence 
and Lu becomes a broad noise component. We only con¬ 
sider the observations for which has good coherence 
(quality factor Q = v^jW > 2). 

For the analysis of the coherent pulsations we use the 
same data selection as for the stochastic analysis. We 
correct the photon arrival times to the Solar System 
barycenter with th e ftool faxbary usin g the optical 
source position of |Hartman et al. (2008), and subse- 
quently correct the da t a tor the binary orbital eph emeris 
(Hartman et al. 2008 2009 Patruno et al. 2012). The 
power spectra oi h'igurel^were made from data corrected 
in this way, and rebinned such that the pulse fundamen¬ 
tal occupies a single frequency bin. In our coherent anal- 
ysis of the pulsations, we fold the data on the pulse period 
( Hartman et al.|2008| ) and for each pulse waveform mea¬ 
sure the sinusoidal amplitudtj^ of the fundamental and 
second harmonic, which we express as a fraction of the 
mean flux. 


3. RESULTS 

In the course of a typical two-week outburst the X- 
ray flux first rises to ~ 70 milliCrab and then falls back 
down. In correlation to this, the upper kHz QPO drifts 
from ~ 300 Hz up to ~ 700 Hz and back down again, and 
so transits the 401 Hz spin frequency twice per outburst. 
We find that at these transits the amplitude of the pul¬ 
sations abruptly changes by a factor ~ 2: it halves when 
in the rise the QPO moves to above the spin frequency 
and then, as illustrated in Figure doubles again on the 
way back. 

We observed three such transitions at high signal-to- 
noise: two where the pulse amplitude approximately 
halves during the outburst rises in 2005 and 2008, and 
one where it app roximately doubles during the decay in 
2002. Figure [^shows the QPO frequency {top frames) 
and the pulse amplitude {bottom frames) across these 
three transitions. In each case the pulse amplitude can 
be seen to systematically drop when the QPO moves to 
frequencies > 401 Hz {blue points)^ and to be system¬ 
atically higher when the QPO is at or below the spin 

^ Note that we report sinusoidal amplitudes, which are a factor 
of \/2 larger than rms amplitudes 



Fig. 1.— Power spectral evolution during the 2002 outburst, 
in chronological order (top-down). The frequency evolution of the 
kilohertz QPO (broad peak) in relation to the spike caused by the 
strictly periodic 401 Hz pulsations is obvious. When the kHz QPO 
frequency (dashed lines) moves through the spin frequency (second 
to third frame) the pulse amplitude (height of the spike) increases 
drastically. Arrows indicate the changes in QPO frequency and the 
initial spike height. 

frequency {red points). In 2002, the observations sam¬ 
ple the source during the actual transition {gray bands). 
The pulse amplitude increases steadily from 3.5 to 5.6% 
on a timescale of hours when simultaneously the QPO 
frequency drifts down from ~ 520 to ~ 400 Hz. Both 
amplitude and frequency change very significantly {17a 
and 8 ct, respectively, see Table [^. In 2005 and 2008 
we observe the opposite transition, with behavior that is 
entirely consistent with that in 2002, but time-reversed. 
The pulse amplitude drops from 7.1 to 4.3% (28cr) and 
6.5 to 3.3% {23a) in 2005 and 2008, respectively, when 
the QPO frequency increases through the 400-520 Hz 
range (I2cr and 5cr change). In all cases, for QPO fre¬ 
quencies < 401 Hz the pulse amplitude remains at the 
high level attained near 401 Hz. 
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2002 2005 2008 




Time (MJD) Time (MJD) Time (MJD) 


Fig. 2.— Evolution of the upper kHz QPO frequency {top frame) and the sinusoidal pulse amplitude {bottom frame) as the QPO 
frequency passes through the 401 Hz spin frequency during outbursts in the years indicated. Observations with the QPO frequency above 
and below the spin frequency are marked in blue and red, respectively. The dashed lines mark the spin-frequency passages and the gray 
bands indicate the full extent of the transits between the two regimes. The double-headed arrows mark pulse amplitude excursions (see 
text). 


TABLE 1 

Spin-frequency Transit Measurements. 



Date 

Exposure 

QPO frequency 

Hi 

H 2 

Outburst 

(MJD) 

(0 

(Hz) 

(%) 

(%) 

2002 

52564.18 

2800 

538 ± 16 

3.49 ± 0.09 

0.76 ±0.09 

52564.97 

3900 

392 ± 10 

5.63 ±0.05 

0.63 ±0.05 

2005 

53524.96 

17500 

397 ±7 

7.10 ±0.07 

1.19 ±0.07 

53525.72 

13900 

518 ±7 

4.35 ± 0.07 

1.35 ±0.07 

2008 

54735.69 

10700 

416 ± 10 

6.49 ± 0.09 

0.98 ±0.09 

54736.44 

12700 

512 ± 18 

3.33 ±0.10 

1.34 ± 0.10 


Note. — Listing {left-to-right)-. the time of observation, total exposure, upper kHz QPO frequency and the fundamental 
(Hi) and second harmonic (H2) sinusoidal pulse amplitudes just before and after each transition for the indicated outbursts. 
The signihcances of these changes were calculated from the errors on the differences using standard error propagation. 


sively 

dHartman et al.||2008 Leahy et al. 2008[ Hartman 

et al.| 

20091 llbragimov & Poutanen||200yi IKajava et al. 

2011 

Morsink & Leahy 2011 Patruno et al. 20121 and 


amplitude. We define an excursion as a sudden pulse 
amplitude increase of 50% or more directly followed by a 
similar amplitude decrease, such that the entire episode 
verifiably takes place within two days. In total there are 


11 such excursions in the data, and we have marked these 
in Figures and P with double-headed gray arrows. 

We find mat the occurrence of these brief excursions 
depends on QPO frequency dichotomously. The excur¬ 
sions only occur when the QPOs are > 401 Hz (see, e.g. 
Figure]^. While they are not limited to any particular 
QPO frequency, the most pronounced excursions cluster 
around the highest QPO frequencies of 600-700 Hz. 

Figure]^ shows the amplitudes of the pulse fundamen- 
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53523 53526 53529 53532 53535 53538 


Time (MJD) 

Fig. 3.— Overview of the 2005 outburst, showing the X-ray flux 
{top frame), the upper kHz QPO frequency (middle frame) and the 
fundamental pulse amplitude (bottom frame) as a function of time. 
Observations with the QPO frequency above and below the spin 
frequency are marked in blue and red, respectively. The dashed 
lines mark the spin-frequency passages and the gray bands show 
the extents of the transits between the two regimes. The double¬ 
headed arrows mark pulse amplitude excursions (see text). The 
other outbursts show similar trends. 

tal {top frame) and second harmonic {bottom frame) ver¬ 
sus upper kHz QPO centroid frequency as we measured 
them in all data across all outbursts. The large step 
in pulse amplitude taking place over a QPO frequency 
range of 400-500 Hz is obvious. When in this transition 
the fundamental amplitude drops from ^ 6 to ~ 3%, 
the weaker second harmonic instead strengthens from 
^ 0.9 to ~ 1.6%. The clustering of the brief excursions 
at high QPO frequencies is also clear among the points 
> 650 Hz. 

The remaining data, at lower signal-to-noise, show be¬ 
havior consistent with that in the three transitions re¬ 
ported above. When in 2005 at MJD 53536.1 the QPO 
moves from above to below 401 Hz, the pulse amplitude 
increases. In 1998 there is no transition; the QPO peak, 
while broad, is always entirely below 401 Hz. Consis¬ 
tent with this, the pulsations are strong and the second 
harmonic is weak. In 2011 kHz QPOs are detected only 
at frequencies > 401 Hz, so no transition is observed. 
In this outburst the pulse amplitude is ^ 4% with, as 
expected, amplitude excursions of 2-3% on top of this. 

4. DISCUSSION 

Our analysis establishes, for the first time, a direct ef¬ 
fect of kHz QPOs on the millisecond pulsations: a large 
change in pulsation properties occurs when the upper 
kHz QPO frequency transits through the spin frequency, 
but reversely a change in pulse amplitude does not pre¬ 
dict the QPO frequency. These findings point toward 



200 300 400 500 600 700 


QPO frequency (Hz) 

Fig. 4. — Fundamental (top frame) and second harmonic (bot¬ 
tom frame) pulse amplitudes as a function of the upper kHz QPO 
frequency of the outbursts as indicated. The dashed line marks the 
spin frequency. 

an origin of the upper kHz QPO in a process involving 
azimuthal motion in the accretion flow that makes the 
accreting plasma interact differently with the neutron- 
star magnetic field depending on whether the azimuthal 
motion is faster or slower than the spin. Although spin 
resonances may play a role (see below), mechanisms rely¬ 
ing on a resonance of the QPO frequency itself with the 
spin cannot explain our observations: the amplitude of 
such a resonance is largest when the QPO frequency is 
near 401 Hz and diminishes away from the resonant fre¬ 
quency. Instead, what is required is a step-wise change 
in the pulsations between QPO frequencies above and 
below the spin frequency. 

Our observation that the pulse amplitude is sensitive 
to whether the QPO frequency is faster or slower than 
the spin immediately suggests an interpretation where 
the QPO is due to orbital motion and the key physics in¬ 
volved is the centrifugal force experienced b y the accret- 


ing plasma as it enters the magnetosphere ( 

Illarionov & 

Sunyaev|1975| 

Ghosh & Lamb| 

19781. In sue 

n a scenario, 
a Keplerian 

plasma outside the magnetosp 

lere orbits in 


disk. Inside the magnetosphere the motions are domi¬ 
nated by the magnetic field and so the plasma is forced 
to corotate with the stellar spin. As long as the radius 
of the magnetosphere, rm, is larger than the corotatiorlj 
radius Tc, the plasma entering the magnetosphere speeds 
up and hence experiences a centrifugal force exceeding 
gravity, inhibiting lateral accretion. For a smaller mag- 

^ The corotation radius r^ is the radius where the Keplerian 

frequency equals the spin frequency: rc = WGM/(2-kVs)\ = 

31 km for an M = 1AMq,Vs = 401 Hz star, where M is the 
neutron-star mass and Ug its spin frequency. 
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netosphere no such centrifugal inhibition of the lateral 
accretion flow will occur and therefore a smaller fraction 
of the accreting plasma is available to flow out of the disk 
plane and toward the poles, causing the pulse amplitude 
to drop. Thus, if the QPO is due to orbital motion at 
the edge of the magnetosphere and hence the QPO fre¬ 
quency can be identified with the Keplerian frequency 
at Tin, a very different accretion pattern is predicted for 
QPO frequencies above and below the spin frequency. 

Another possible mechanism that can explain our re¬ 
sult can apply to any model in which the QPO is due 
to an azimuthal motion at the inner edge of the disk. 
This motion could be Keplerian, but might also be, e.g., 
preces sional, or that of a n azimuthally propagating disk 
wave ( van der Klis|2006 and references therein). As the 
inner edge of the disk rnoves inward it will pass through 
a resonance radius, r^, where the QPO frequency equals 
the spin frequency and a resonance may occur. For even 
smaller inner disk radii the QPO frequency increases fur¬ 
ther, but the resonance may continue to exist in the disk 
at r^. It is at least logically possible that such a res¬ 
onance present in the disk flow at radius leads to a 
different accretion pattern at radii inside r^. How pre¬ 
cisely the resonance affects the pulse amplitude remains 
unspecihed in this scenario, which in this generic sense 
could formally even work for an azimuthally symmetric 
QPO mechanism, and/or by interaction with the pulsar 
beam instead of the magnetic field. 

In the above two scenarios, the accretion pattern 
changes when the inner edge of the disk passes through 
or r^, respectively. A third possibility is that the relevant 
radius is ri„. It has been suggested that it is possible for 


the Keplerian disk to penetrate the magnetosphere and 
continue its orbital motion inside Cm, with t he inner edge 


of th e disk at rin set by radiative stresses (Miller et al. 


19981, so that r-m can be either larger or smaller than 


r^. If the QPO originates from the inner edge of the 
disk, and its frequency is near the spin frequency as ri„ 
passes through r^, differences in magnetic coupling as¬ 
sociated with Tin being larger or smaller than r„i might 
conceivably explain the change in the pulsations. For a 
QPO caused by Keplerian orbital motion these condi¬ 
tions might be fulfilled, as during outb urst SAX J1808 


may be near spin-equilibrium (r,„ ~ rc, Hartman et al. 


20081 . 

So) although some more complex explanations cannot 
formally be excluded, the interpretations discussed here 
predominantly point to a QPO that is produced by az¬ 
imuthal motion at the inner edge of the accretion disk, 
most likely orbital motion, and a change in the interac¬ 
tion of the accreting plasma with the neutron star mag¬ 
netic field near a QPO frequency equal to the spin fre¬ 
quency. The most straightforward scenario involves cen¬ 
trifugal inhibition of the accretion flow that sets in when 
the upper kHz QPO becomes slower than the spin. This 
scenario implies that the QPO is produced by orbital 
motion at the edge of the magnetosphere as it passes 
through 401 Hz. 


We would like to thank the anonymous referees whose 
comments improved this work. This research has made 
use of data obtained through the High Energy Astro¬ 
physics Science Archive Research Center Online Service, 
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